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ABSTRACT 


A  study  was  conducted  to  determine  the  effect  of 
specimen  size  on  the  mechanical  response  of  rock.  Specimens 
of  Cedar  City  tonalite  (Cedar  City,  Utah)  and  Charcoal  Black 
granite  (Cold  Springs,  Minnesota)  ranging  in  size  from  2  in, 
dia.  to  32  in.  dia.  (22  in.  dia,  for  the  tonalite)  were  tested 
in  triaxial  compression.  In  addition  to  the  standard  triaxial 
test  stress  trajectory,  one  specimen  of  each  rock  type  was 
loaded  under  uniaxial  strain  to  simulate  true  one -dimensional 
compression.  Test  data  included  axial  and  circumferential 
strain  at  up  to  30  locations  on  the  largest  specimens,  surface 
wires  to  indicate  crack  propagation  on  the  12  inch  and  larger 
specimens,  axial  and  radial  stresses  and  the  confining  pressure 
media  temperature.  Data  were  recorded  under  loading,  unloading 
and  reloading  conditions.  The  test  data  for  each  specimen 
size  and  rock  type  were  fit  to  a  non-linear  hysteretic  model 
to  determine  variation  in  the  model  parameters  as  a  function 
of  specimen  size. 
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TRIAXIAL  TESTS  ON  LARGE  ROCK  SPECIMENS 


1.  INTRODUCTION 

Currently  a  number  of  structures  are  being  located 

1  2 

underground,  many  of  which  are  in  rock.  It  is  anticipated  * 
that  the  number  of  such  underground  structures  will  increase 
enormously  in  the  next  two  decades.  It  may  be  concluded  that 
the  number  of  structures  in  rock  will  increase  correspondingly. 
In  spite  of  this  anticipated  underground  construction  activity, 
little  is  known  about  the  mass  behavior  of  rock.  It  is 
qualitatively  recognized,  of  course,  that  the  rock  mass  behaves 
distinctly  differently  from  the  matrix  in  small  samples,  but 
this  effect  of  size  is  not  well  understood,  and  cannot  be 
predicted,  to  date,  with  any  degree  of  confidence.  These 
variations  are  attributed  to  joints  and  other  weakness  planes 
in  the  rock  mass.  Yet,  most  of  the  test  procedures  for 
determining  strength  and  elastic  moduli  information  about 
rocks  are  based  on  these  properties  on  small  samples3"6. 

It  is  both  difficult  and  expensive  to  try  to  perform  properties 
tests  on  larger  specimens.  Hence  some  relationships  between 
the  various  sizes  need  to  be  established.  It  was  the  intent 
of  this  program,  to  study  the  mechanical  behavior  of  rock 
specimens  of  different  sizes,  and  to  determine  how  these 
relate  to  each  other. 

During  this  program,  tw:>  varieties  of  rock,  Charcoal 
Black  granite  from  Cold  Springs,  Minnesota  and  Cedar  City 
tonalite  from  Cedar  City,  Utah,  were  investigated.  Four 
sizes  of  rock  cores,  2-in.  dia.,  4-in.  dia.  12-in.  dia., 
and  32-in.  dia.  (22-in.  in  case  of  the  tonalite)  were  subjected 
to  triaxial  compression  tests.  The  specimens  were  strain 
gaged  to  determine  the  elastic  moduli!.  One  test  on  the 
large  core  samples  of  each  rock  type  was  uniaxial  strain 
only.  The  data  obtained  was  correlated  to  ascertain  trends. 
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Accumulation  of  sufficient  data  on  size  effects 
with  various  rock  types  may  be  expected  to  finally  lead 
to  predictive  equations  for  the  rock  mass.  It  should  be 
borne  in  mind,  however,  that  the  strength  properties  of  smaller 
rock  samples  are  governed  by  microflaws  such  as  crystal 
lattice  disturbances,  grain  boundaries,  minute  voids  (pores) 
and  weak  cementing  material.  The  strength  of  the  rock  mass 
is  dependent  on  macroflaws,  i.e.  failure  planes,  joints, 
cavities,  and  intrusions.  It  also  appears  probable  that 
the  rock  mass  is  more  affected  by  environmental  conditions 
than  the  matrix.  In  general,  environmental  factors  induce 
a  time-dependent  reduction  in  strength  by  chemical  alteration, 
rheological  deformation,  and  structural  failure.  The  roles 
of  residual  stresses  in  the  rock  matrix  and  orogenic  forces 
on  the  mass  are  not  well  understood;  hence  their  effect  on 
rock  behavior  requires  further  study.  The  presence  of  water 
significantly  affects  the  behavior  of  the  rock  mass  as  well 
as  the  matrix.  In  general  it  tends  to  reduce  the  structural 
strength,  but  whether  the  manner  in  which  this  is  accomplished 
in  the  two  cases  is  the  same,  has  not  been  clearly  established. 
Finally,  it  should  be  mentioned  that  the  configuration  of 
the  rock  mass,  and  the  base  material  on  which  it  rests  probably 
controls  its  gross  behavior.  It  would  be  difficult  to  simulate 
these  conditions  in  laboratory  samples,  but  the  overall  effects 
maybe  estimated  by  varying  size  and  end  conditions. 
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2.  PREVIOUS  STUDIES 


2.1  Effect  of  Size 


Several  studies  have  been  conducted  to  determine  the 
effect  of  size  of  specimens.  Most  of  these  have  been  with 
unconfined  compressive  tests.  Rice  and  Enzian^  showed  that 
the  compressive  strength  of  2  1/2  -  4  in.  coal  cubes  was  about 
2500  psi,  whereas  a  54  in.  cube  of  the  same  coal  failed  at 

Q 

300  psi.  Greenwald,  Howard,  and  Hartmann  related  the  size 
of  coal  pillars  to  their  compressive  strength  at  an  experimental 
face  by  over-cutting  and  side-cutting.  There  was  a  decrease 
in  strength  with  increasing  pillar  height  relative  to  the 
cross-sectional  area,  the  breaking  stress  o  for  small  pillars 
approximately  satisfying  the  relationship 


a 


c 


©c 


0.5 


» 


where  b  is  the  lateral  dimension  of  the  pillar  and  h  is  its 
height.  Burton  and  Phillips9  studied  the  relationship  between 
size  and  compressive  strength  of  cubes  of  anthracite  over  the 
range  1/4  In.  to  4  In.  and  found  that 


o  «  a -0.45*3. 05 

r.  * 


where  a  is  the  linear  dimension  of  the  cubes.  Millard,  Newman 
and  Phillips  extended  the  range  of  these  observa tions ,  and 
used  Griffith's  theory  of  crack  propagation  to  explain  the 
failure  phenomenon.  Both  these  investigations  by  Phillips  and 
his  associates  may  be  said  to  show  that  approximately 

ac°C  a"°‘5» 

Evans  and  Pomeroy11  found  that  the  relation  between  the  mean 
crushing  load  (<p)  and  the  side  of  the  cube  (a)  is  'of  the 
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form 


<P  <* 


where  fl  ■  constant  varies  between  1.5  and  2,  i.e., 


Holland12 
follows : 


a 


oc  a 


-0.5  to  0 


c 

has  suggested  a  coal  pillar 


a 


design  formula  as 


where  a  -  least  width  of  pillar  (inches),  h  -  thickness  of 
pillars  (inches) ,  and  k  -  coefficient  depending  on  type 
of  coal.  This  is  valid  for  a/h  ratios  between  1  and  10. 

For  rocks,  the  types  of  relationships  developed  are 
slightly  different.  According  to  the  Weibull  theory1  , 


0€ 


(Zs) 


where  o  -  tensile  or  compressive  strength  of  the  rock  from  a 
standard  laboratory  test, 
om  -  equivalent  strength  of  the  rock  mass, 

v  ■  volume  of  the  rock  mass, 

v  -  volume  of  the  test  sample,  and 

04  -  constant  (with  values  near  10  for  rocks). 

The  relation  established  by  Protodyakonov1^  was  of  the  type 


1  + 


-  u 

s  +  a 


where  s  -  spacing  between  major  discontinuities  in  the  rock 
mass,  e.g.,  joints,  beds. 

a  -  dimension  of  the  test  specimen,  usually  diameter 
for  compression  tests,  and 
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8  -  mass  fracture  coefficient  (in  compression: 1-2  for 
igneous  rocks,  1-3  for  competent  sedimentary 
rocks,  3-10  for  weak  rocks;  in  tension:  approximately 
double  these  values) . 

Grobbelaar^  based  on  the  work  of  Epstein^,  Bieniawski^, 
and  others,  found  that  the  formulae  relating  the  modal 
strength  of  the  weakest  elements  and  its  standard  deviation, 
based  on  the  weakest  link  theory,  are 

°N  "  °u  ‘  °s  <21°8n)0,5-  1/2  [log(logN)  +  log (4  )}  (21ogN)-° 

and  o  (N)  -  c  xir(121ogN)'0,5 
s  s 

or  os(N)  -  os(N0)  (logN0/logN)°-5 

where  N  -  number  of  flaws  in  the  large  cubic  specimen, 

NQ  -  number  of  flaws  in  the  small  cubic  specimen  or 
unit  cube, 

aN  -  modal  strength  of  the  weakest  link  in  a  sample 
containing  N  elements, 

ou  -  average  strength  of  a  unit  cube  of  material 
containing  Nu  elements , 

o8  -  standard  deviation  of  the  modal  strength  of 
samples  containing  NQ  elements , 
and  afl(N)  ■  standard  deviation  of  the  modal  strength 
of  the  weakest  element  in  a  sample  containing 
N  elements. 

These  formulae  are  based  on  the  "weakest  link  theory", 
which  can  be  analyzed  mathematically  if  it  is  assumed  that 
the  frequency  of  occurrence  of  events  is  a  continuous 
function  (e.g.  Weibull^  or  normal  distribution).  The  "links" 
in  this  case  are  the  macroflaws  (or  cracks)  in  the  bulk 
material;  not  the  microflaws. 
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Glucklich  and  Cohen^®*^  have  indicated  that  effects 
other  than  statistical  exist,  since  the  total  stored  elastic 
energy  increases  with  specimen  volume.  The  energy  released 
at  onset  of  fracture  is  related  to  initiate  fracture;  in  other 
words,  this  reduces  the  rock  strength.  This  phenomenon  has 
been  recently  discussed  by  Baecher 

2.2  Effect  of  Confinement 

There  have  been  numerous  studies  investigating  the 

effects  of  various  aspects  of  confinement  on  rocks.  It  is 

not  intended  to  review  all  of  these  completely  in  this  section. 

Most  of  the  pertinent  work  has  been  briefly  discussed  by 
21 

Swanson  .  The  earliest  experimental  work  was  performed  by 
22  23 

Adams  and  von  Karman  .  However,  significant  headway  was 

not  made  until  the  initiation  of  work  by  Griggs^  and 
25-26 

his  coworkers  D  .  Since  then  a  number  of 

researchers  have  conducted  various  types  of  studies  under 

pressure  several  of  which  were  presented  at  a  symposium  on 

rock  deformation37.  Baidyuk3®  has  summarized  some  of  the 

Russian  and  American  work.  Research  in  this  area  is  still 
21  29  30 

very  active  *  *  .  All  of  this  work  has  been  performed 

with  small  rock  specimens,  a  few  inches  in  diameter.  As  a 

result  considerable  light  has  been  shed  on  the  behavior  of 

the  rock  matrix  and  the  criteria  of  failure.  Refinements 

to  the  Griffith  hypothesis  have  been  proposed3*”33  and 

appear  to  explain  the  rock  fracture  process  under  confinement 

fairly  well.  The  extrapolation  of  these  theories  to  larger 

rock  masses  is  of  doubtful  value  and  hence  large  scale 

field  testing  has  to  be  resorted  to  .  The  U.  S.  Bureau  of 

Mines  has  undertaken  a  rather  comprehensive  program  to 

collect  field  data  with  the  intention  of  correlating  it  into 

a  hypothesis35.  The  contributions  of  Hoek36,  Bieniawski37, 

7  38  39 

Cook  ,  Wawersik  and  HoupertJ7  to  the  mechanism  of  brittle 
failure  in  rock  deserve  to  be  noted  even  though  the  studies 
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were  not  conducted  under  a  confined  state  of  stress. 

To  the  best  of  the  authors'  knowledge  no  experimental 
investigations  on  large  rock  samples  under  confinement 
have  been  performed. 
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3.  EXPERIMENTAL  APPARATUS 


In  order  to  conduct  triaxial  tests  for  the  range  of 
specimen  sizes  used  on  this  project,  four  triaxial  cells 
were  set  up  as  shown  below: 


Chamber  Specimen  Maximum 

I.D.  (in.)  Diameter  (in.)  Chamber  Pressure 

(ksi) 


Maximum 
Axial  Load 
lb.  x  106 


4.0 

2.06 

6.5 

3.65 

14.7 

10-12 

48.3 

22-32 

30  0.375 

30  0.990 

20  3.40 

20  axial  36 . 5 
10  confining 


In  a  standard  triaxial  cell  the  axial  load  is  supplied  by  an 
external  loading  machine.  However,  in  order  to  achieve  the 
large  end  loads  required  for  the  tests  in  this  program,  these 
chambers  were  separated  into  two  regions  by  sliding  pistons. 
The  general  configuration  is  shown  in  Fig.  1.  One  region 
contained  the  specimen,  and  was  pressurized  to  the  desired 
confining  pressure.  Axial  load  was  transmitted  to  the  rock 
by  the  sliding  piston.  The  maximum  axial  stress  in  the  rock 
depends  upon  the  ratio  of  the  rock  and  piston  areas  and  the 
difference  in  the  confining  pressure  and  the  axial  chamber 
pressure.  The  specimen  stresses  are  given  by  the  following 
equations: 

°2  "  °3  "  P3 
AOl  -  <P1*P3>  Jr 

where 

Ap  ■  piston  area 
Ar  ■  specimen  area 
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Chamber  Closure  with 
Instrumentation  Lead- 
Out 


Confining  Pressure 
Region 


Test  Specimen 


Sliding  Piston  with 
Loading  Blocks 


Axial  Pressure  Region 


Fig.  1  TYPICAL  TRIAXIAL  CELL  SCHEMATIC 
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1 L  -  axial  chamber  pressure 
"confining  chamber  pressure 
AO.  -  deviator  stress 

*  i 

al*  a2  an<*  a3  "  Pr*-nc*-Pal  stresses. 

i  :  :  I 

! 

Reference  to  Figure  1  and  the  above  equations  confirms  that 
if  P^  -  P^,  A aim  0  and  the  specimen  is  uhdet  hydrostatic 
stress  (a^  -  ■  0^)  .  For  P^  ■  0,  the  specimen  is  unconfined 

with  (?2  -  0^  m  0  and  ■  P^  Ap/Ar.  '■  1 

3.1  Small  Test  Cells  ,  '  1 

:  •'  ! 

The  three  smaller  test  cells  are  incorporated  into 
a  testing  system  with  centralized  controls,  instrumentation, 
and  pumping  systems.  A  schematic  of  this  system  is  shown 
in  Fig.  2.  The  tests  conducted  in  this  system  consisted  of 
initial  hydrostatic  loading  up  to  the  desired  confinihg 
pressure,  followed  by  triaxial  compression  at  constant  ,c3 
above  that  pressure.  At  least  one  load-unload-reload  cycle 
was  observed  for  each  test.  Provision  was  made  for  pressure 
cros 8 -conned tions  between  the  confining  pressure  and  axial 
pressure  chamber  volumes  to  insure  hydrostatic  conditions 
during  the  hydrostatic  test  phases.  , 

,  1  1 

3 . 2  48-Inch  Test  Cell  ( 

I 

This  test  cell  is  shown  in  Figs.  3  and  4.  The  basic 
unit  is  a  48”  I.D.  by  86”  working  length , chamber  having  1 
20,000  psi  design  working  pressure.  The  chamber  walls'  are 
built  up  from  rings  12"  long  which  are  held  in  place  by  a 
3/4"  thick  liner  bn  the  inside  diameter.  The  entire  axial 
load  is  carried  by  a  flexible  reaction  frame  which  was  built 
up  from  steel  strap.  The  l4o  ton  weight  of  thia  chamber  is 
substantially  less  than  the  weight  that  would  have  been 
required  by  conventional  chamber  design.  This  is  the  1 

i 

1  !  '  > 
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SCHEMATIC  OF  SYSTEM  FOR  2  IN.,  4  IN.,  AND  12  IN.  TRIAXIAL 
SPECIMENS 


NOT  REPRODUCIBLE 


Fig. 


48  IN.  DIA.  CELL  WITH  32  IN.  DIA.  SPECIMEN  BEING 
LOADED 


largest  chamber  currently  available  at  IIT  Research  Institute, 
and  is  capable  of  applying  axial  loads  of  36  million  pounds 
to  rock  specimens  as  large  as  3  1/2  ft.  in  diameter. 

As  can  be  seen  in  the  schematic,  the  pumping 
and  control  systems  for  this  unit  are  simpler  than  for 
the  small  chambers.  A  separate  pump  was  used  at  each  end 
of  the  chamber.  Accumulators  were  not  used  because  the 
chamber  volume  itself  is  large  in  comparison  with  other 
available  pressure  chambers. 

The  operation  of  this  chamber  is  similar  to  that 
of  the  smaller  chambers,  except  that  the  turn-around  time 
between  tests  in  on  the  order  of  a  week  instead  of  an  hour. 

The  tests  in  this  chamber  differed  slightly  from  those  in 
the  smaller  test  cells.  In  order  to  maintain  seal  integrity, 
a  positive  pressure  differential  of  about  200  psi  was  maintained 
across  the  sliding  piston  during  the  '’hydrostatic"  portions 
of  the  triaxial  tests,  and  the  axial  pressure  was  not  allowed 
to  drop  below  400  psi  at  the  bottom  of  the  load-unload-reload 
cycle.  Had  these  precautions  not  been  taken,  there  would 
have  been  danger  of  upsetting  the  piston  seal,  thus  aborting 
the  remainder  of  the  test. 

Ttoo  of  the  large  specimens  were  loaded  to  produce 
uniaxial  strain.  Figure  5  shows  the  two  transducers  were 
designed  to  provide  measurements  of  average  radial  strain. 

One  transducer  consisted  of  a  full  strain  gage  bridge  using 
manganin  wire  elements.  The  active  elements  each  consisted 
of  4  turns  of  manganin  wire  wound  tightly  around,  but  not 
bonded  to,  the  specimen.  Temperature  and  pressure  compensation 
were  achieved  by  similar  elements  loosely  wound  on  the 
specimen  next  to  the  active  elements.  The  second  transducer 
consisted  of  an  open  steel  loop  pinned  to  the  specimen. 

The  open  ends  of  the  loop  were  jointed  by  a  flexible  strap 
mounting  foil  strain  gages  as  shown  in  Figure  5b.  Small 
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-  Active  Elements 
Passive  Elements 


5a  -  Manganin  Wire  Gages 


Fig'  5  STRAIN  GAGES  FOR  UNIAXIAL  STRAIN 

iilibTb 


deflections  at  the  diametrically  opposed  pins  cause  relatively 
large  bending  strains  in  the  flexible  strap.  This  transducer 
suffered  from  somewhat  erratic  operation,  as  is  usual  with 
any  point-to-point  gages  on  rock  specimens. 

3.3  Specimen  Preparation 

Specimens  of  the  Charcoal  Black  Granite  were  bought 
from  Cold  Spring  Granite  Co.  in  Cold  Spring,  Minnesota,  but 
the  Cedar  City  tonalite  was  supplied  by  DASA.  In  both  cases, 
the  larger  cores  (10  or  12-in.  and  22  or  32-in.  dia.)  were 
obtained  as  such  from  the  source,  but  the  2-in.  and  4-in.  dia. 
specimens  were  cored  in  the  laboratory  from  extra  rock  obtained 
in  each  case.  These  smaller  cores  were  cut  parallel  to  the 
axis  of  the  larger  cores  so  as  not  to  introduce  complications 
because  of  anisotropy. 

End  preparation  of  the  2-in.  and  4-in.  dia.  cores 
consisted  of  facing  on  a  lathe  until  the  ends  were  plane  and 
parallel  to  0.001  Inches.  The  larger  cores  were  capped  in 
a  specimen  cage  to  permit  handling.  The  capping  material 
used  was  a  steel-filled  epoxy.  Figures  3  and  7  each  show 
an  assembled  32  inch  specimen  in  its  cage.  The  cage  tie 
rods  were  designed  with  end  fittings  that  would  accept 
tensile  load  only.  Since  the  maximum  tensile  load  that 
could  be  applied  to  the  cage  by  these  tie  rods  corresponded 
to  30  psl  compressive  stress  in  the  rock  specimen,  the  effect 
of  the  cage  on  the  rock  was  negligible. 

An  array  of  foil  strain  gages  were  mounted  on  each 
specimen.  The  number  used  ranged  from  three  rosettes  on 
the  2  in.  cores  to  thirty  rosettes  on  the  32  in.  cores. 

These  were  two  element  rosettes  with  1/4  in.  gage  length 
placed  with  the  direction  of  rolling  parallel  to  the  specimen 
axis.  The  gage  placement  procedure  included  the  following  steps 

*  grind  the  rock  surface 

*  apply  two  thin  coats  of  gage  cement 
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*  visual  inspection  for  voids  in  the  cement  base 

*  affix  and  wire  gage 

*  apply  two  coats  of  gage  cement  for  water  proofing 

*  check  gage  for  continuity  and  response  (soft 
eraser)  and  replace  if  necessary. 

In  several  cases  it  was  necessary  to  move  a  rosette  slightly 
away  from  its  pre -determined  location  because  of  local  flaws 
or  drill  shot  embedded  in  the  surface  (especially  in  the  case 
of  the  tonalite) . 

In  addition  to  the  gages,  break  wires  were  cemented 
to  the  surface  of  the  12-in.  and  greater  diameter  specimens. 

This  wire  was  No.  36  manganin  wire  well  bonded  to  the  surface 
so  that  the  wire  would  break  if  a  crack  propagated  across 
it  during  a  test.  Figure  6  shows  typical  instrumentation  for 
a  12 -in.  core. 

The  instrumented  cores  were  waterproofed  with  latex 
cement  over  the  foil  gages,  and  at  least  two  coats  of  latex 
paint  with  a  thickening  agent  spread  over  the  entire  rock  surface. 
Waterproofing  is  very  important  in  a  triaxial  test  since 
both  the  loading  conditions  and  the  character  of  the  rock 
can  be  changed  by  intrusion  of  oil  into  the  rock  voids.  In 
those  tests  where  the  specimens  could  not  be  loaded  to 
failure  on  the  second  load  cycle,  the  specimen  was  recovered 
intact  and  stripped  of  paint  for  visual  inspection.  In 
several  cases,  traces  of  oil  were  found  under  the  paint,  but 
typically  the  specimens  were  completely  dry.  In  no  case  was 
there  enough  oil  to  do  more  than  dampen  a  very  small  area 
on  the  specimen  surface. 

3.4  Instrumentation  and  Data  Reduction 

The  instrumentation  on  this  program  included  foil 
8 train  gages,  temperature  transducers,  pressure  gages  and 
break  wires  on  the  larger  specimens.  All  instrumentation 
except  the  pressure  gages  were  recorded  using  an  automatic 
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data  acquisition  system.  The  pressures  were  read  using 
Hiese  Super  Accurate  bourdon  gages  and  inserted  manually 
onto  the  data  system  output  as  they  were  read. 

Due  to  the  bulk  of  data  involved,  the  data  were 
reduced  and  plotted  using  the  1108  Univac  computer  at  IITRI. 
In  addition  to  printed  output,  the  following  plots  were 
produced  for  each  test: 

Shear  strain  vs.  deviator  stress 
Volumetric  strain  vs.  mean  stress 
Poisson's  ratio  vs.  mean  stress 
Octahedral  shear  stress  vs.  mean  stress 
Circumferential  and  axial  strain  vs.  axial 
stress 

Elastic,  shear  and  bulk  moduli  vs.  mean  stress. 

Data  from  each  test  was  then  assembled  by  hand  to  show  trends 
from  test  to  test. 
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4.  EXPERIMENTAL  PROGRAM 

The  program  Included  triaxial  tests  on  specimens 
ranging  from  2-in.  to  32-in.  dia.  One  uniaxial  stiain  test 
was  conducted  on  the  largest  available  specimen  for  each 
rock  type,  these  being  a  32 -in.  dia.  granite  and  a  22 -in. 
dia.  tonalite.  The  uniaxial  tests  were  conducted  by  monitoring 
mean  radial  strain  and  maintaining  sufficient  confining  pressure 
to  hold  this  strain  as  near  zero  as  possible.  The  experiment::  1 
program  is  summarized  in  Table  1.  1ITRI  was  able  to  procure 
the  granite  core  from  Cold  Springs  Quarry,  Minnesota.  The 
tonalite  specimens  were  delivered  by  DASA  which  was  unable 
to  obtain  more  than  one  large  specimen  of  tonalite. 

The  loading  and  unloading  sequence  in  the  large 
chamber  is  shown  by  the  photographs  in  Figs.  7,  8  and  9. 

Typical  failures  are  seen  in  Figs.  10  thru  13. 
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TEST  NUMBERS  IN  EXPERIMENTAL  PROGRAM 


I 


I 


i 


i  i  ■ 

I 

,  NOT  REPRODUCIBLE 


Fig.  7  LOADING  SPECIMEN  INTO  48^IN,  DIA. 
TRIAXIAL  CELL 


I 


I 


NOT  REPRODUCIBLE 


Fig.  8 


INSERTING  SLIDING 
TRIAXIAL  CELL 


PISTON  INTO  48-IN. 


DIA. 
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not  reproducible 


Fig.  9  REMOVING  FRACTURED  SPECIMEN  FROM  48 -IN. 
DIA.  TRIAXIAL  CELL 


NOT  REPRODUCIBLE 


Fig.  10  32-IN.  DIA.  CHARCOAL  BLACK  GRANITE  SPECIMEN 

AFTER  TESTS  27-28  (a  -  5000  psi) 

c 


26 


NOT  REPRODUCIBLE 


Fig.  11 


32— IN.  DIA.  CHARCOAL 

AFTER  TESTS  29-30  (a 

'  c 


BLACK  GRANITE  SPECIMEN 
■  10,000  psi) 


NOT  REPRODUCIBLE 


Fig.  12  32-IN.  DIA.  CHARCOAL  BLACK  GRANITE  SPECIMEN 

AFTER  TEST  31  (UNIAXIAL  STRAIN) 


NOT  REPRODUCIBLE 


Fig.  13  22-IN.  DIA.  CEDAR  CITY  TONALITE  SPECIMEN 

AFTER  TEST  32  (UNIAXIAL  STRAIN) 
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5 .  DISCUSSION  OF  RESULTS 


•5,1  Theoretical  Considerations  for  Constitutive  Equations 

The  large  amount  of  data  generated  during  this  program 
makes  a  simple  report  of  experimental  data  inappropriate. 

The  reduced  data  Itself  consists  of  an  8  1/2  inch  thick  stack 
of  printed  computer  output.  Since  it  was  obvious  at  the 
start  of  the  program  that  this  would  be  the  case,  parallel 
computer  output  in  the  form  of  graphical  plots  was  provided 
for  inclusion  in  the  report  as  an  appendix.  In  order  to 
summarize  the  data  for  discussion  in  the  text  of  the  report, 

It  was  felt  that  the  most  convenient  form  would  be  a  comparison 
of  the  test  data  with  a  model.  The  model  parameters  can 
then  be  evaluated  as  a  function  of  specimen  size  much  more 
easily  than  the  actual  data. 

5.1.1  Formulation  of  Constitutive  Equations 

Concepts  from  the  theory  of  plasticity  have  recently 
been  used  to  describe  the  behavior  of  particulate  materials 
such  as  soil  and  rock.  If  a  stress  space  is  constructed  with 
coordinate  axes  representing  the  three  principal  stresses, 
tne  state  of  stress  of  an  element  of  material  is  described  by 
a  point  in  this  stress  space.  The  locus  of  all  possible  stress 
states  which  cause  yielding  of  the  material  is  a  surface  called 
the  yield  or  limit  surface.  This  surface  separates  the 
stress  states  attainable  in  the  material  from  those  which  are 
not  attainable.  An  elastic-plastic  material  is  one  which 
behaves  as  an  elastic  material  at  stress  states  within  the 
limit  surface,  and  as  a  plastic  material  at  stress  states  on 
the  limit  surface.  A  material  model  thus  must  include  both 
a  description  of  the  limit  surface  and  the  elastic  equations 
that  hold  within  the  limit  surface. 
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5.1.2  Limit  Surface 


The  limit  surface  In  principal  stress  space  is  usually 
expressed  by 

Vj|~  -  f(Jx)  a) 

In  which  Is  the  first  invariant  of  tbi  stress  tensor  and 
j£  is  the  second  invariant  of  the  deviatoric  stress  tensor. 
These  invariants  are  defined  by 


-  ax  +  o2  +  a3 
(a1-c3'>2  +  (ol  - 


°2>2  +  <a2  ‘ 


(2) 

(3) 


in  which  o^,  cr2 j  and  a 3  are  the  three  principal  stresses.  An 
equivalent  method  of  describing  the  yield  surface  is  to  write 
the  octahedral  shear  stress,  r  t,  as  a  function  of  the  octahedral 
normal  stress,  a  ,  or 

GC  w 


Toct  "  f<aoct) 


(4) 


in  which 


and 


°oct  -  1  (01  +  52  +  °3> 


roct  "  f 


(a^Og)2  +  (a^"a3)2  +  (°2~°3)2 


(5) 


1  1/2 


(6) 


One  form  of  a  yield  surface  commonly  used  is  the 
Mi8e8'  yield  criteria  given  by 

J2  ‘  Ko  (7) 

in  which  Kq  is  the  yield  limit  of  the  material  in  simple  shear. 
Equation  (7)  describes  the  surface  of  a  cylinder  in  principal 
stress  space  with  the  axis  inclined  equally  to  all  three 
principal  stress  axes  as  shown  in  Fig.  14.  According  to  Mises' 
yield  criterion  the  shear  strength  of  a  material  is  independent 
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Fig.  14  YIELD  SURFACE  FOR  MISES'  YIELD  CRITERION 


(8) 


of  the  mean  stress  a  ,  where 

m 


In  granular  materials,  Including  most  soils  and  rocks,  It 
has  been  established  that  shear  strength  Increases  with 
mean  stress  due  to  frictional  behavior.  For  these  materials 
the  Coulomb-Mohr  yield  surface  Is  more  applicable,  and  Is  given 
by 

Vi]"  -  A  +  B  J1  (9) 

In  which  A  and  B  are  constants.  Equation  (9)  describes  the 
surface  of  a  cone  In  principal  stress  space  with  the  same 
axis  as  the  Mlses'  yield  criteria  shown  In  Fig.  14. 

5.1.3  Elastic  Equations 

Constitutive  equations  for  an  Isotropic  linear  elastic 
material  are  given  by 

alj  “  *6ijekk  +  2^eij  (10> 

In  which 

ay  ■  stress  tensor 
■  8 train  tensor 
Cjj  -  Kronecker's  delta 
are  Lame's  constants. 

It  is  sometimes  convenient  to  separate  the  volumscric  response 
which  is  produced  by  mean  stress  from  the  distortions!  response 
which  is  produced  by  shear  stress.  This  is  accomplished  through 
the  deviator  strain  tensor,  e'jj*  and  the  devia top  stress 
tensor,  which  are  defined  as 

f  ij  "  eij  ~  1  ekk  6ij  (W 

°ij  "  °ij  *  J  °kk  sij  (12) 
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ekk  +  2G‘ij 


(13) 


The  constitutive  equations  (IQ)  may  then  be  rewritten  in  the 

•  :■  1  i 

form 

“ij  ■ k  ®ij 

in  which  k  and  G  are  Respectively  the  bulk  modulus  and  shear 
modulus.  These  moduli  are  defined  by  ■  , 

a,  , 

tk 

3(ej+£2+£3) 


k  "  7 


°-7 


kk 

ekk 

£ 

£  ij 


(14) 


(15) 


Nonlinear  behavior  may  be  incorpora jted  by  using  an  incremental 
constitutive  relation  and  defining  k  and  G  as  functions  of 
one  or  more  of  the  stress  Invariants.  Inelastic  behavior 
may  be  Incorporated  by  using  separate  valuer  of  k  and  G  for 
loading  and  unloading.  The  unloading  moduli  must  be  greater 
than  the  loading  moduli  to  prevent  energy  generating  hystersis 
loops .  i 

t  1 

5.2  Typical  Behavior  of  Rock 

5.2.1  Hydrostatic  Behavior 


Past  work  has  shown  that  the  stress -strain  behavior 
.  •  1 
of  rock  subjected  to  hydrostatlq  pressure  is  non-linear.  A 

typical  curve  of  mean  stress  versus  volumetric  strain  from 

Paulding*®  is  shown  in  Figure  15.  In  an  unloaded  state  a 

rock  specimen  usually  has  an  initial  porosity,  including  some 

small  cracks.  As  hydrostatic  pressure  increases,  the  cracks 

begin  to  close  and  the  porosity  decreases,  causing  the  rock 

to  become  stiffer  and  less  compressible.  At  higher  pressures 

the  response  becomes  linear  and  compressibility  of  the  rock 

may  be  predicted  from  the  individual  compressibilities  of 

its  constituent  minerals.  The  Initial  porosity  may  be  determined 

by  extending  the  linear  portion  of  the  curve  back  to  the 

strain  axis.  Based  on  data  presented  by  Clark41,  shown  in 

Fig.  16,  it  appears  that  the  influence  of  the  iriitial 
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Hydrostatic  Pressure 


f 


I 


J 


Compressibility  of 
Solid  Material 


Fig.  16  COMPRESSIBILITY  OF  GRANITE,  GABBROS,  AND 
DIABASES:  FROM  CLARK41 


porosity  on  compressibility  becomes  negligible  at  a  pressure 
of  about  2000  bars  or  29  ksi. 


The  behavior  of  granular  materials  under  hydrostatic 
pressure  has  been  studied  analytically  through  the  use  of 
models  consisting  of  various  packings  of  elastic  spheres^. 
Hertz's  contact  theory^  has  been  used  and  has  shown  that 
the  volumetric  strain  is  directly  proportional  to  mean  stress 
to  the  two- thirds  power.  This  nonlinearity  in  response  is  due 
to  the  geometry  of  the  packing  of  the  individual  particles 
and  not  the  properties  of  the  material  comprising  the  spheres. 
The  Hertz  theory  may  be  applicable  if  rock  may  be  considered 
as  a  granular  material. 


A  nonlinear  hysteretic  model  used  by  Seaman  and 

A  A 

Whitman  to  study  the  behavior  of  sand  appears  to  be  suitable 
for  representing  the  hydrostatic  behavior  of  rock.  The 
stress -strain  curve  for  this  model  is  shown  in  Fig.  17.  For 
virgin  loading 


a 


m 


(16) 


and  for  unloading  and  reloading 

°m  ‘  A2  <ev  -  M>"  <17> 

in  which 

a  «  mean  stress 
m 

ey  *  volumetric  strain 

e»g  ■  residual  volumetric  strain 
A-^  and  A2  ■  material  properties 
This  model  represents  both  the  nonlinear  and  the  inelastic 
behavior  of  the  rock  observed  under  hydrostatic  loading. 

Application  of  the  mathematical  model  in  a 
computer  code  would  perhaps  be  most  convenient  in  an 
incremental  form  using  tangent  values  of  bulk  modulus, 
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kt.  This  modulus  is  a  function  of  mean  stress  and  for  a 
given  stress  represents  the  slope  of  the  curve  at  that 
stress  level. 


For  loading 

.  do 
k.  ■  _ m 

v 


„l/n  „n-l/n 
■  n  a,  O" 
l  m 


and,  for  unloading  and  reloading 

kt  •  £»  - » 4/n  C1/n 


(18) 


(19) 


Note  that  in  this  model  the  modulus  for  unload-reload  it, 
larger  than  the  modulus  for  virgin  loading  at  any  given 
stress  level. 

5.2.2  Triaxial  Compression 

A  triaxial  compression  test  normally  consists  of  two 
phases.  First,  a  hydrostatic  confining  stress,  a  ,  is  applied 
to  the  specimen  so  that  the  principal  stresses  are  all  equal- 
to  oQ.  Then,  two  principal  stresses,  a2  and  o3,  are  kept 
constant  at  oq  while  the  third  principal  stress,  o^,  is 
increased.  For  these  conditions  one  of  the  components  of  the 
de via tor  stress  is 


all  "  °i  "  1/3  (al  +  2a3>  "  2/3  (°i“a3>  (20) 

If  the  material  is  isotropic  two  principal  strains,  e? 

and  e3,  are  also  equal,  meaning  that  the  corresponding  component 

of  deviator  strain  is 


£{l  "  ei  “  l/3  (£i  +  2e3)  -  2/3  (e^^  -  e3)  (21) 

Thus,  for  triaxial  conditions  the  shear  modulus  may  be 
determined  from 


qU  °1  ‘  g3 

~  Eg) 
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Idealized  behavior  of  a  rock  specimen  loaded  ii< 

21 

triaxial  compression  as  presented  by  Swanson  is  shown 
in  Fig.  18.  The  relation  between  deviator  stress  and  strain 
is  linear  up  to  yield.  At  stress  states  beyond  yielding  the 
relation  becomes  nonlinear  and  plastic  behavior  occurs. 

The  modulus  for  unloading  is  equal  to  that  for  virgin  loading 
to  the  yield  point. 

Volume  change  behavior  is  the  same  as  that  for 
hydrostatic  loading  at  stress  states  below  yield.  At  yield 
the  volume  of  the  specimen  begins  to  increase  even  though 
the  mean  stress  is  compressive  and  also  increasing.  This 
tendency  for  rock  to  expand  when  yielding  is  termed  dilatancy 
and  is  characteristic  of  almost  all  particulate  materials. 

At  low  stress  levels  shear  modulus  increases  relatively 
rapidly  with  mean  stress;  however,  at  higher  stresses  it 
remains  relatively  constant.  Torsional  wave  velocities 
were  measured  in  cylinders  of  different  granites  exposed 
to  hydrostatic  pressure  between  1  kg/cm^  and  4,000  kg/cm^ 
by  Birch  and  Bancroft*1^.  From  the  wave  velocities  the  modulus 
of  rigidity,  which  is  identical  to  the  shear  modulus,  was 
determined  and  is  shown  as  a  function  of  mean  stress  in  Fig.  19. 
Also  shown  in  Fig.  19  are  values  of  shear  moduli  for  specimens 
of  Charcoal  Black  granite  and  Cedar  City  tonalite  determined 
from  triaxial  tests  performed  during  this  study.  The  moduli 
determined  by  Birch  and  Bancroft^  and  those  for  the  Charcoal 
Black  granite  are  in  good  agreement  and  increase  approximately 
with  the  1/10  power  of  the  mean  stress.  Values  of  shear 
modulus  for  the  Cedar  City  tonalite  are,  however,  radically 
different.  At  mean  stress  levels  below  about  5  ksi,  the 
shear  modulus  increases  with  the  5/6  power  of  mean  stress. 

At  higher  stress  levels,  the  shear  modulus  increases  less 
rapidly.  One  probable  cause  for  this  change  in  behavior 
is  the  relatively  high  initial  porosity  of  the  tonalite. 
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Fig.  18 
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Values  of  shear  modulus  shown  in  Fig.  19  are  tangent 
values,  Gt,  which  refer  to  the  slope  of  the  deviator  stress - 
strain  curve  at  a  given  stress  vtlue.  For  nonlinear  behavior 
the  tangent  shear  modulus  is  a  function  of  mean  stress  and 
may  be  calculated  from 


1/2 


~  o3) 

“  £3) 


(23) 


The  functional  relationship  with  mean  stress  may  be  written 

<»> 

t  m 

in  which  c  is  a  constant. 

5.2.3  Uniaxial  Strain 

The  uniaxial  strain  test  is  a  special  type  of  triaxial 
test  in  which  the  lateral  strains,  62  and  e^,  are  equal  to 
zero  throughout  the  test.  One  method  of  obtaining  this 
condition  is  through  monitoring  the  lateral  strains  and 
controlling  the  lateral  stresses,  Oj  and  o^,  to  maintain 
zero  strain.  In  soil  mechanics,  for  the  uniaxial  strain 
loading,  the  ratio  between  the  lateral  and  axial  stress  is 
called  the  coefficient  of  earth  pressure  at  rest,  KQ,  where 

K  -  ^3  .  (25) 

°1 

Data  from  this  test  may  be  used  to  compute  bulk  and  shear 
moduli  from  the  following  equations : 

G  ■  1/2  -=■  (1  -  K  ) ,  and  (26) 

fcl 

k  =  1/3  ^  (1  +  2K0).  (27) 

5.2.4  Poisson's  Ratio  and  Young's  Modulus 

Two  constants  are  necessary  to  describe  the  stress- 
strain  behavior  of  an  isotropic  elastic  material.  One  set 
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of  constants,  k  and  G,  which  separate  volumetric  and 
deviatoric  behavior  have  been  described  previously. 

Another  pair  of  constants  which  are  commonly  used  are 
Poisson's  ratio,  v,  and  Young's  modulus,  E.  The  constants 
E  and  v  are  derived  directly  from  a  trlaxial  test  in  which 
the  lateral  stresses,  and  0^,  are  equal  to  zero.  For 
this  special  test  condition. 

E  -  i ,  and  (28) 

el 

v  -  -  (29) 

el 


For  more  general  test  conditions 


E 


9kG 


3F+T? 


and 


(30) 


3k  -  2G 
v  -  2(31*4)  * 


(31) 


The  use  of  tangent  values  of  bulk  and  shear  moduli,  kt  and 
Gt,  in  equations  (30)  and  (31)  will  result  in  tangent  values 
for  Poisson's  ratio,  vt,  and  Young's  modulus,  Et. 


For  hydrostatic  loading  the  deviatoric  stresses  and 
strains  are  zero  and  G  is  undetermined.  Therefore.  E  and  v 


are  also  indeterminate  from  a 
test  with  a ^  equal  to  zero,  E 
equations  (28)  and  (29).  For 
unequal  to  zero  but  constant, 
may  be  computed  from 

d  a* 

Ec  *  aq'  > and 


hydrostat.  For  a  trlaxial 
and  v  may  be  computed  from 
a  trlaxial  test  with  0^ 
tangent  values  of  E  and  v 

(32) 

(33) 
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For  more  general  conditions  equations  (30)  and  (31)  must 
be  used. 


For  uniaxial  strain  tests  the  following  expressions 
may  be  used  .  .  K  ) 

— n-Tin — -• and  (») 


0, 

E  «  _JL 


v-  rnr 


(35) 


5.2.5  Yield  and  Fracture  Criteria 

In  principal  stress  space  the  restraint  02  ■  03 
defines  a  plane  which  includes  the  o^-axis  and  the  axis 
of  the  yield  surface  shown  in  Fig.  14.  As  02  -  ©3  is  a 
condition  maintained  during  the  triaxlal  test  all  stress 
states  in  this  test  lie  in  the  ©2  *  ©3  plane.  Thus,  yield 
and  fracture  criteria  may  be  described  by  a  curve  in  this 
plane.  Curves  representative  of  different  failure  criteria 
are  shown  in  Fig.  20  with  plotted  vs.  Y  on  orth°gonal 
axes.  The  stress  path  for  a  constant  03  triaxial  test  is 
also  shown  in  Fig.  13. 

A  thorough  discussion  of  yield  and  fracture  surfaces 
as  they  apply  to  the  behavior  of  rock  has  been  given  by 
Swanson^.  The  yield  surface  defines  all  stress  states 
at  which  the  rock  ceases  to  be  an  elastic  material.  At 
stress  states  within  the  yield  surface  some  form  of  elasric 
equations  govern  the  material  behavior.  At  states  on  the 
yield  surface  equations  from  plasticity  govern  the  behavior. 
It  has  been  shown ^  that  most  rocks  exhibit  strain -hardening 
behavior;  i.e.,  the  increase  of  stress  beyond  initial  yield 
causes  the  yield  surface  to  move  outward.  Continued  increase 
of  stress  causing  outward  movement  of  the  yield  surface 
will  eventually  result  in  failure  of  the  rock.  The  locus 

of  all  1  tress  states  causing  failure  is  called  the  fracture 
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Fig.  20  YIELD  CRITERIA 


I 


I 


i  surface  or  failure  surface.  , 

.  *  /  1  !  .' 

5.3,  Analysis  of  Experimental  Results 

Experimental  results  are  presented  in  a  form 
consistent  with  the  elastic-plastic  constitutive  model 

I  ; 

presented  earlier.  Elastic  properties  are  presented  as 
functions, of  mean  stress  and  are  given  in  terms  of  tarigent 
values  which  are  compatible  with  a  nonlinear ; incremental 
form  of  elasticity  equations.  Yield  and  fracture  criteria 

j  i 

are  presented  in  the  form  Ji  ■=  ( J, ) . 

,  ’  1  1  L  i  1  ' 

5.3.1  Charcoal  Black  Granite 

I  , 

Typical  volume  change  data  from  a  cbnstant  confining 
i  pressure  triaacial  test  is  shown  in  Fig.  21  along  with  a  curve, 
representing  the  nonlinear  hy^teretic  model  fitted  to  the  > 

,  ■  1  i 

experimental  data.  The  constants  for  the  mbdel  were  conven¬ 
iently  determined  from  logarithmic,  plots  of  the  data,  In1 
the  middle  portion  of  the  curves  experimental  results  and 

.  I  !  j 

the  model  are  in  close  agreement.  At  relatively  low  stress 
levels,  less  than  2  ksi,  the  model  prfedicts,  stresses  ttyat 
are  smaller  than  those  measured.  Fracture  of  this  specimen 
occurred  at  a  mean  stress  of  17  ksi  and  yielding  had  begun 
at  about  14  ksi  which  results  in  dilatant  behavior  of  the 
specimen  at  these!  stress  levels,  'fhus,  as  would  be  expected 
the  model  and  experimental  results  dp  differ  at  stress  levels 
greater  than  14  ksi  for  this  particular  specimen. 

I  !  1  ' 

Tangent  values  of  bulk  modulus,  Kt,  were  determined 
from  the  experimental  data  and  plotted  as  a  function  of 
mean  stress  in  Figs.  22  and  23.  These  values  represent 
the  computed  slope  between  actual  consecutive  data  points 
rather  than  the  slbpe  of  a  smooth  curve  drawn  through  the 
data  points.  This  procedure  introduces  "noise"  i,r«  the 
plots  since  the  data  points  are  rather  close  together  and 
the  Uncertainity  in  the  computed  volumetric  strain  is  rather 

large  in  comparison  with  the  change  in  volumetric  strain  over 
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Mean  Stress 


Bulk  Modulus,  k,_,  ksi  Bulk  Modulus,  k*.,  ksi 


Bulk  Modulus,  k.  ,  ksi  Bulk  Modulus,  k  ,  ksi 


khim  umn 


Fig.  23  BULK  MODULUS  VS.  MEAN  STRESS  -  12”  0  AND  32”  0  CBG 

50 


the  short  interval  between  data  points.  Each  plot  is 
superimposed  on  the  values  of  bulk  modulus  predicted  by 
the  nonlinear  hysteretic  model  for  the  particular  core 
size.  Several  trends  are  immediately  apparent.  First, 
the  experimental  data  fits  the  model  best  in  the  hydrostatic 
range  of  each  test,  and  departs  from  the  model  at  low 
stress  levels  and  at  high  stress  levels  where  yielding  is 
to  be  expected.  Second,  the  fit  to  the  model  is  best  for 
the  small  specimens  and  worse  for  the  large  specimens. 

This  is  probably  a  result  of  the  relatively  early  slippage 
that  took  place  in  the  larger  specimens.  The  larger  specimens 
began  slippage  at  low  levels  of  deviator  stress,  so  that 
earlier  departure  from  the  model  is  to  be  expected.  Of  course, 
it  is  to  be  remembered  that  the  model  was  developed  for  small 
rock  specimens. 

Table  2  is  a  summary  of  the  model  parameters  for  the 
tests  on  Charcoal  Black  granite.  The  smaller  cores  behaved 
quite  consistently,  and  the  model  parameters  were  simply 
averaged  to  produce  the  final  parameters  for  each  test. 

The  larger  tests  were  less  consistent,  particularly  with 
respect  to  the  values  for  A^.  The  final  parameters  for  the 
32-in.  specimens  were  obtained  from  a  composite  log-log 
plot  of  data  for  tests  27,  28  and  30. 

An  interesting  comparison  may  be  made  between  the 
analysis  of  an  ideal  packing  of  elastic  spheres  as  discussed 
earlier  and  the  observed  hydrostatic  behavior  of  the  rock 
specimens.  Values  of  Mn"  in  the  nonlinear  hysteretic  model 
vary  from  1.59  to  1.93  with  an  average  of  1.74  for  the 
specimens  of  Charcoal  Black  granite.  The  corresponding 
values  of  V  for  the  ideal  packing  is  1.5.  For  granular 
soils  n  values  between  1.5  and  2.0  have  been  observed^. 
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NONLINEAR  HYSTERETIC  MODEL  PARAMETERS  FOR  TEST  ON  CHARCOAL  BLACK 

GRANITE 
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Typical  deviatoric  stress  and  strain  data  from  a 
constant  confining  pressure  triaxial  test,  «  5  ksi, 
are  shown  in  Fig.  24.  The  curves  are  very  nearly  linear 
at  stresses  less  than  approximately  30  ksi  which  is  the 
level  at  which  initial  yielding  of  the  specimen  occurs. 

Also  the  unload  curve  is  nearly  elastic  with  a  very  small 
residual  strain.  Similar  data  was  obtained  from  tests  at 
other  confining  pressures  except  the  unconfined  compression 
test,  a ^  ■  0.  For  this  test  the  slope  of  the  deviatoric 
stress  vs.  deviatoric  strain  curve  increased  slightly  as  the 
stress  increased  meaning  that  the  shear  modulus,  G,  also 
increases  with  stress  at  low  stress  levels. 


Tangent  values  of  shear  modulus,  Gfc,  were  determined 
from  the  experimental  data  and  are  shown  in  Fig.  24  (b) , 

25  (b),  26  (b)  and  27  (b)  for  specimen  diameters  of  2  in., 

4  in.,  12  in.  and  32  in.  respectively.  The  stress  paths 
of  the  particular  tests  considered  for  each  diameter  are 
shown  in  Figs .  25  (a) ,  26  (a) ,  27  (a) ,  and  28  (a) .  As 
observed  from  the  data  the  shear  modulus  increases  with  mean 
stress  especially  at  the  low  stress  levels.  An  equation 
of  the  form 


Gt  ■  co 
t  m 


1/10 


(36) 


in  which  c  is  a  constant,  was  fitted  to  the  experimental 
results  and  is  shown  in  Fig.  25  through  28.  As  relatively 
little  inelastic  behavior  was  observed  in  the  shear 
deformations  the  model  for  shear  modulus  is  the  same  for 
both  loading  and  unload-reload.  Model  predictions  are  in 
agreement  with  observed  behavior  except  at  low  stresses. 


A  summary  of  the  mathematical  models  for  bulk 
and  shear  moduli  along  with  the  constants  determined  from 
the  experimental  results  is  presented  in  Table  3.  Values 
of  Young's  modulus  were  computed  from  equation  (30)  using 
the  predicted  values  of  bulk  and  shear  modulus  and  are  shown 
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Fig.  24  DEVIATOR  STRESS  VS.  STRAIN,  CHARCOAL  BLACK  GRANITE 
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Test  12 


Test  14,  34 


b)  Shear  Modulus  and  Young's  Modulus 

Fig.  25  TANGENT  MODULI  FOR  CHARCOAL  BLACK  GRANITE,  2 -IN.  DIA. 
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b)  Shear  Modulus  and  Tangent  Modulus 


Fig.  27  TANGENT  MODULI  FOR  CHARCOAL  BLACK  GRANITE,  12-IN.  DIA. 
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8  12 
Mean  Stress 
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Table  3 


II  UW* 


ELASTIC  PROPERTIES  OF  CHARCOAL  BLACK  GRANITE 


Bulk  Modulus 

loading,  kt  -  nAx1/n  am  ^ 
unload  -  reload,  kfc  -  nAg1^  ly^n 

Shear  Modulus 

loading,  Gt  ■  C  Gfc  and  om  in  ksi 

unload  -  reload,  Gfc  ■  C  am^^»  Gt  and  o>m  in  ksi 


Young's  Modulus 


E 


9  kt  G 
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as  a  function  of  mean  stress  in  Fig.  25  through  28.  Poisson's 
ratio  values  were  computed  from  equation  (31)  and  are  shown 
in  Figs.  29,  30,  and  31  for  the  different  specimen  diameters. 
Again,  except  at  stresses  below  1  ksi  the  model  fits  the  exper¬ 
imental  data. 

5.3.2  Cedar  City  Tonallte 

The  nonlinear  hysteretic  model  was  fitted  to  the 
various  core  sizes  of  Cedar  City  tonalite  by  use  of  log-log 
plots  of  mean  stress  versus  volumetric  strain.  The  resulting 
model  parameters  were  used  to  compute  bulk  moduli  in  the 
same  fashion  as  was  done  for  the  Charcoal  Black  granite. 

Figures  32,  33  and  34  show  the  actual  bulk  moduli  as  a 
function  of  mean  stress.  The  solid  symbols  indicate  data 
from  the  hydrostatic  portions  of  the  various  tests,  the  open 
symbols  the  triaxial  compression  portions.  In  general,  the 
hydrostatic  data  fall  near  the  model  and  the  triaxial  data 
depart  from  the  model.  The  data  points  do  not  scatter 
symmetrically  about  the  model  because  the  model  was  fit 
to  the  volumetric  strains  rather  than  the  bulk  moduli. 

The  same  observations  apply  to  this  data  as  to  the 
corresponding  Charcoal  Black  granite  data.  The  model  for 
each  core  size  fits  best  in  the  intermediate  stress  ranges, 
and  the  data  scatter  increases  with  specimen  size.  The 
nonuniformity  of  this  rock  type  is  particularly  noticeable 
in  the  12-in.  cores.  Tests  17  and  38  had  quite  different 
moduli  below  5  ksi,  but  appeared  to  converge  to  the  same 
stiffness  above  7  ksi.  This  may  be  expected  if  the  existing 
cracks  are  closing  in  this  pressure  range. 

The  scale  effect  is  shown  in  Figs.  35  and  36.  Figure 
35  is  a  plot  of  the  mean  stress  versus  volumetric  strain 
predicted  by  the  model  for  the  three  specimen  sizes.  Figure  36 
is  a  common  plot  of  the  bulk  moduli.  The  major  difference  between 
the  2-in.  and  4-in.  models  is  that  the  4  in.  model  is  slightly 
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Fig.  31  POISSON'S  RATIO  VS.  MEAN  STRESS  FOR  32"  0  CHARCOAL 
BLACK  GRANITE 
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Fig.  32  BULK  MODULUS  FOR  2"  CEDAR  CITY  TONALITE 
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Fig.  33  BULK  MODULUS  FOR  4"  CEDAR  CITY  TONALITE 
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Fig.  35  VOLUMETRIC  STRAIN  FOR  CEDAR  CITY  TONALITE  MODELS 
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a  Lillet  below  3  ksi.  This  results  in  a  horzontal  displacement 
of  the  curves  in  Fig.  35.  If  Fig.  35  were  extended  to  higher 
stress  levels,  the  2-in.  and  4-in.  curves  would  cross  as  a 
result  of  the  higher  2-in.  bulk  modulus.  The  actual  differences 
in  these  curves  are  probably  a  result  of  specimen  variation 
rather  than  scale  effects.  The  12-in.  specimens  appear  to  be 
significantly  less  stiff  than  the  smaller  specimens. 

Figure  37  shows  the  scale  effect  upon  shear  modulus 
Gt  and  elastic  modulus  Et.  Data  from  all  tests  are  included, 
and  again  the  scatter  is  caused  by  the  short  interval  between 
data  points.  It  can  be  seen  that  the  2-in.  and  4-in.  specimens 
were  not  significantly  different,  and  that  the  12-in.  specimens 
were  less  stiff  at  lower  mean  stress  values.  The  shear  modulus 
did  not  fit  the  one-tenth  power  of  mean  stress  model,  so  fair 
lines  were  passed  through  the  data  by  eye.  The  elastic 
modulus  model  was  computed  using  values  of  the  bulk  modulus 
models  in  Figs.  32,  33  and  34,  and  appears  to  represent 
the  experimental  elastic  moduli  well. 

Figure  38  shows  plots  at  Poisson's  ratio  for  the 
various  size  specimens.  The  typical  behavior  is  to  start 
out  somewhat  below  the  model  and  rapidly  increase  to  a 
value  greater  than  would  be  predicted  by  this  model.  The 
tests  having  15  ksi  confining  pressure  maintain  a  value 
of  approximately  125%  of  the  model  value  for  a  period  before 
yielding  begins,  while  the  tests  at  lower  confining  pressures 
appear  to  maintain  a  mono tonic  increase.  As  the  yield  point 
approaches,  the  value  of  Poisson's  ratio  may  increase  because 
of  the  dilatancy  of  the  specimen. 

5.3.3  Uniaxial  Strain  Tests 

Uniaxial  strain  tests  were  conducted  on  two  specimens, 
a  22-in.  dia.  Cedar  City  tonalite  specimen  and  a  32-in.  dia. 
Charcoal  Black  granite  specimen.  The  results  from  this  type 
of  test  should  be  applicable  to  unidirectional  compression 
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Fig.  37  TANGENT  MODULI  FOR  CEDAR  CITY  TONALITE 
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Fig.  38  POISSON'S  RATIO  FOR  CEDAR  CITY  TONALITE 


of  an  element  in  a  rigidly  confined  mass  or  plane-wave 
propagation.  The  radial  strains  were  monitored  during  the 
test  by  the  auxiliary  averaging  strain  gages  described 
in  the  experimental  apparatus  section  of  this  report. 
Confining  pressure  was  controlled  to  maintain  these  gage 
readings  as  near  zero  as  possible. 

Figure  39  shows  the  axial  and  lateral  strains  and 
the  stress  trajectories  for  both  tests.  As  can  be  seen, 
the  radial  strains  did  go  tensile  by  about  250  microstrains 
in  test  31  and  100  microstrains  in  test  32.  The  Cedar  City 
tonalite  (test  32)  failed  at  an  axial  stress  of  20  ksi. 

The  behavior  of  the  Charcoal  Black  granite  was 
somewhat  more  complex,  primarily  because  this  specimen  was 
able  to  support  load  after  slippage,  and  the  test  was 
continued  for  two  full  cycles.  The  offset  on  the  unload 
cycle  at  11,5  ksi  was  caused  by  a  gage  failure.  The  offset 
on  the  reload  cycle  at  26.5  ksi  was  accompanied  by  a  major 
rock  shift. 


Figures  40  and  41  show  the  shear  strain  and  volumetric, 
strain  for  both  tests.  Note  that  for  this  test  condition 
there  appears  to  be  little  difference  between  the  two  rock 

£  £  j 

types.  Both  have  shear  moduli  Gt  between  4  x  10  and  5.x  10  „ 
Any  curvature  in  the  Charcoal  Black  granite  shear  modulus 
is  swallowed  by  the  initial  zero  offset,  which  may  or  may 
not  be  real.  The  bulk  modulus  for  the  Cedar  City  tonalite 
increases  from  1.6  x  10^  to  5.2  x  10^  ksi  as  mean  stress 
goes  from  zero  to  9  ksi.  The  Charcoal  Black  granite  bulk 
modulus  ranged  from  1.7  x  10  ksi  at  1,5  ksi  to  4  x  10°  ksi 
at  a  mean  stress  of  1.0  ksi.  Corresponding  values  for  elastic 
modulus  and  Poisson’s  ratio  are  listed  below:  ' 
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5.4  Strength  Properties 

Figure  42  indicates  the  maximum  stress  levels  attained 
by  each  specimen.  Due  to  a  less  favorable  intensification  ratio 
in  -he  14-in.  triaxlal  cel! ,  there  were  no  failures  of  the  12 -in. 
specimens.  Future  programs  would  need  to  employ  higher  stress 
levels  to  insure  failure. 

The  32-in.  dia.  Charcoal  Black  granite  specimens  reached 
maximum  stress  levels  of  about  one  half  that  of  the  2-ln.  specimens. 
The  actual  scale  effect  is  probably  somewhat  greater.  Slippage 
along  a  joint  plane  will  produce  instant  failure  in  a  2-in. 
specimen,  while  a  32-in.  specimen  not  only  survives  minor  slippage, 
but  can  be  loaded  to  greater  stress  levels. 

The  Cedar  City  tonalite  was  much  more  variable.  As 
might  be  expected  from  a  weathered  rock,  the  strength  at  low 
mean  stress  levels  is  quite  low.  This  corresponds  to  a  low  "c" 
value  on  a  plot  of  the  Mohr-Coulomb  failure  envelope. 

Enough  data  was  not  obtained  in  this  program  to  apply 
any  statistical  theories,  such  as  those  of  Weibull13  or 
Grobbelaar^3.  It  is  hoped  this  will  become  possible  when  more 
data  is  accumulated. 
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Fig.  42  MAXIMUM  VALUE  OF  V^T"  REACHED  FOR  ALL  TESTS 


Over  the  course  of  this  program,  27  triaxial  tests 
were  run  on  samples  ranging  from  2  in.  to  32  in.  dia.  This 
is  the  first  time  that  specimens  of  this  size  have  ever 
been  tested  under  laboratory  conditions.  Because  of  the 
large  volume  of  data  collected,  the  results  of  this  program 
are  presented  in  term  of  a  nonlinear  hysteretic  model  proposed 
by  Seaman  and  Whitman^  and  having  the  form 


°m  “  Alev  loadinS 

°m  "  A2(Vevi)n  unloadln8 


where 


a  is  mean  stress 
m 

ev  is  volumetric  strain 

e...  is  residual  volumetric  strain. 


The  data  fit  this  model  well  at  intermediate  stresses. 
However,  the  model  predicts  that  bulk  modulus  is  zero  when 
mean  stress  equals  zero,  and  bulk  modulus  goes  to  infinity 
at  infinite  mean  stress.  In  fact,  rocks  display  a  finite 
bulk  modulus  at  zero  mean  stress  and  a  constant  bulk  modulus 
under  high  stresses.  In  spite  of  this,  the  data  depart 
from  the  model  only  below  stresses  of  about  3  ksi  and 
at  higher  stresses  when  yielding  occurs. 

The  model  parameters  as  fitted  to  each  core  size 
are  summarized  in  Tables  4  and  5.  Up  to  12-in.  dia.  it 
seems  the  influence  of  specimen  size  on  the  mechanical 
properties  is  small.  There  is  not  enough  data  on  the  larger 
samples  to  draw  definite  conclusions.  The  largest  scale 
effect  was  found  in  the  Cedar  City  tonalite  at  mean  stress 
levels  below  5  ksi.  Both  the  model  and  the  experimental 
data  indicated  low  bulk  modulus  and  elastic  modulus  in  this 
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range,  although  these  parameters  Increase  rapidly  with  mean 
stress. 

The  Influence  of  specimen  size  on  strength  Is 
pronounced  on  both  rock  types.  The  32 -In.  Charcoal  Black 
granite  specimens  failed  at  stress  levels  about  one  half 
those  of  the  2-in.  specimens.  The  mode  of  failure  also 
was  markedly  different.  The  2-ln.  specimens  failed 
catastrophically,  while  the  32-in.  specimens  failed  Initially 
by  a  series  of  small  slips. 

With  the  accumulation  of  more  test  data  It  may 
become  possible  to  develop  a  model  which  more  clearly 
predicts  the  behavior  of  the  large  rock  cores.  It  may 
also  become  meaningful  to  verify  statistical  hypotheses. 

It  should  be  noted  that  with  the  larger  cores,  the 
rock  depicted  a  load  carrying  capacity  even  after  failure 
had  Initiated.  This  is  often  observed  In  the  field.  Besides, 
the  mechanism  of  failure  was  different  In  these  larger 
cores.  This  is  significant  in  that  at  about  30-in.  dla. 
or  larger  perhaps  the  type  of  rock  failure  obtained  appears 
to  be  similar  to  that  observed  in  the  field.  Perhaps  models 
developed  for  this  size  can  predict  results  closer  to  field 
observe  tlooA. 
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APPENDIX  A  -  COMPUTER  PLOTS  OF  TEST  DATA 

Data  was  reduced  and  plotted  using  the  Univac  1108  at 
1ITRI.  This  appendix  includes  the  three  most  significant  plots 
for  each  test,  these  being: 

shear  strain  vs.  deviator  stress 

volumetric  strain  vs.  mean  stress 

axial  and  lateral  strains  vs.  axial  stress. 

In  several  cases,  e.g.  test  19  a  vs.  e„,  the  data  went  off- 
scale  during  the  plotting  subroutines.  When  this  occurred,  the 
data  was  plotted  at  an  arbitrary  value  within  the  plot  scale. 

The  plots  are  arranged  in  order  of  rock  type,  core 
diameter,  confining  stress  and  test  number. 
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